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We describe how sideband cooling techniques, prevalent in quantum optics, may be applied to 
large spin ensembles in magnetic resonance. Using the Tavis-Cummings model in the presence of 
a Rabi drive, we solve a Markovian master equation describing the joint spin-cavity dynamics to 
derive cooling rates as a function of ensemble size. Our calculations indicate that a spin ensemble 
containing roughly 10^^ electron spins may be polarized to a non-thermal equilibrium state in a 
time many orders of magnitude shorter than the typical thermal relaxation time. The described 
techniques permit the efficient removal of entropy for spin-based quantum information processors 
and fast polarization of spin samples. The proposed application of a standard technique in quantum 
optics to magnetic resonance also serves to reinforce the connection between the two fields, which 
has only recently begun to be explored in detail due to the development of hybrid designs for 
manufacturing noise-resilient quantum devices. 



I. INTRODUCTION 

Efficient removal of entropy from a quantum system 
presents a significant challenge toward the development 
of new quantum technologies and devices. High purity 
quantum states that may be quickly initialized and re- 
set are necessary for the application of quantum error 
correcting codes to suppress and mitigate the effects of 
noise and errors that naturally occur in quantum infor- 
mation processors, sensors, and communication devices 
[ij . Additionally, in spectroscopic applications the signal- 
to-noise ratio increases significantly with state purity, al- 
lowing for the detection of small spin ensembles. 

A spin ensemble may be naively prepared in a pure 
state by simply moving to low temperatures, where ther- 
mal excitations are not energetic enough to cause sig- 
nificant transitions out of the ground state. However, 
the required temperatures are often impractical to ob- 
tain or require sophisticated and expensive equipment. 
Additionally, the time required for the spin system to 
reach thermal equilibrium with the environment - the 
energy relaxation time, Ti - often becomes very long at 
low temperatures, limiting the rate at which spin resets 
and signal averaging may be applied [2]. 

A variety of techniques for removing entropy from a 
quantum system are commonly used, including dynamic 
nuclear polarization (DNP) [2|l3|, algorithmic cooling [4], 
optical pumping [5,, laser cooling |6]-t8j, and microwave 
cooling pun]) among others. Recently, it was demon- 
strated that superconducting qubits may be prepared in 
an arbitrary pure state through sideband cooling by a 
high quality factor (high-Q) cavity [TTl [T^] . We discuss 
in this work how similar microwave cooling techniques 
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should also be applicable to ensemble spin systems in 
magnetic resonance, despite the relatively small coupling 
between the cavity and a single spin. In particular, we 
present a theoretical model for how a high-Q resonator 
(cavity) may be used to actively drive an ensemble spin 
system to a highly pure, non-thermal equilibrium state 
on a timescale that is significantly shorter than the ther- 
mal Ti. 

II. MATHEMATICAL MODEL 

We consider an inductively driven ensemble of non- 
interacting spin-1/2 particles quantized in a large static 
magnetic field and magnetically coupled to a high-Q cav- 
ity. In the presence of the drive the spins interact with 
the cavity via coherent radiative processes and may be 
treated quantum mechanically as a single collective mag- 
netic dipole coupled to the cavity [T^. In analogy to 
quantum optics, we describe the spin-cavity dynamics as 
being generated by the Tavis-Cummings (TC) Hamilto- 
nian [Ml [15] . Assuming the control field to be on reso- 
nance with the Larmor frequency of the spins, the spin- 
cavity Hamiltonian under the rotating-wave approxima- 
tion (RWA) is given by H = Hq + H^it) + Hj, with 

= cJcO,^ a + oJsJz (2-1) 
Hj^{t) = nRCOs{ujst)J^ (2.2) 
Hi = g{aKj^+aJ+), (2.3) 

where a^(a) are the creation (annihilation) operators de- 
scribing the cavity, Q.^ is the strength of the drive field 
(Rabi frequency), ujc is the resonant frequency of the cav- 
ity, LOg is the Larmor resonance frequency the spins, and 
g is the coupling strength of the cavity to a single spin 
in the ensemble in units oi h = 1. Here we have used 
the notation that J a = J2f=i^a^ the total an- 

gular momentum spin operators for an ensemble of Ng 
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spins. The TC Hamiltonian is also known as the Dicke 
model [16! and has been studied extensively for quantum 
optics (for a recent review see [T7]). 

The eigenstates of Hq are the tensor product of photon- 
number states for the cavity and spin states of collec- 
tive angular momentum in the Jz direction: J, "niz) 
Here n = 0, 1, 2, . . ., = —J, — J -I- 1, . . . , J — 1, J, 
and J = Ns/2. The collective excitation number of 
the joint system is given by Nex — a^a + (J^ -I- J). 
The interaction term Hj commutes with N^x, and 
hence preserves the total excitation number of the sys- 
tem. It drives transitions between the state |n)^|J, m^)^ 
and states \n + l)JJ,mz — 1) ^ and |n — 1)^| J, -t- 1)^ 
at a rate of \/{n + 1)( J( J + 1) — mz{mz — 1) and 
y^n{J{J + 1) — mz{mz + 1), respectively. 

After moving into an interaction frame defined by 
Hi ~ ijjs{a^a + Jz), the spin-cavity Hamiltonian is trans- 
formed to 



^(1) 



Hi 



(2.4) 



= Suja^a + ^Jx+ g{a^J- + aJ+) (2.5) 



for atomic ensembles, but not restricted to the single- 
excitation manifold |21j - enhances spin polarization at 
a rate that may exceed the thermal relaxation rate. 

We note that the spin-cavity exchange coupling also ex- 
ists in the absence of the Rabi drive, and permits cooling 
of the spin system by matching the resonance frequency 
of the spin system to the cavity resonance. However, 
this process is thermally driven, and thus corresponds to 
incoherent radiative processes that do not display the co- 
herent collective enhancement obtained in the presence 
of a drive [T3] . This Purcell effect in magnetic resonance 
systems has been previously noted and is normally small 
enough to be neglected due to the low probability of spon- 
taneous emission in spin systems [22l |23j . 

To model the cavity-induced cooling of the spin system 
we use an open quantum system description of the cavity 
and spin ensemble. The joint spin-cavity dynamics may 
be modelled using the time-convolutionless (TCL) master 
equation formalism [24j, allowing the derivation of an 
effective dissipator acting on the spin ensemble alone. 

The evolution of the spin-cavity system is described by 
the Lindblad master equation 



where 5uj — ujc — uJs is the detuning of the drive from 
the cavity resonance frequency and we have made a sec- 
ond RWA to remove the doubly rotating terms in the 
interaction frame. If we now move into a second inter- 
action frame of H2 = Suia^a + flj^ Jx/^, the Hamiltonian 
transforms to 



H+nJt) 

'J X 



(2.6) 



¥)*atj(-) 



Hon At) = g[e-''^t)a + e''^'A) 
H^nAt) = '{{ 

(x) 

where Jj_ = Jj, ± iJz are the spin-ladder operators in 
the a;-basis. 

In analogy to Hartmann-Hahn matching in magnetic 
resonance cross-relaxation experiments [T8H20] for Suj > 
we may set the cavity detuning to be close to the Rabi 
frequency of the drive, so that A = Soj — r2i?/2 is small 
compared to Suj. Then, by discarding the high frequency 



terms in Eqn (2.6) the interaction Hamiltonian reduces 
to the H-n.R flip-flop exchange interaction between the 
cavity and spins in the x-basis: 



H 



(2.7) 



From here we will drop the (x) superscript and just note 
that we are working in the Jx eigenbasis. 

Isolating the spin-cavity exchange interaction allows 
efficient energy transfer between the two systems, per- 
mitting them to relax to a joint equilibrium state in the 
interaction frame of the control field. The coherent en- 
hancement of the ensemble spin-cavity coupling - sim- 
ilar to the enhancement of the vacuum Rabi frequency 



dt 



p{t)^Ci{t)p{t)+V,p{t) 



(2.8) 



where £/ is the super operator Ci{t)p — —i[Hi{t), p] 
describing evolution under the interaction Hamiltonian 
( |2.7[ ) , and Vc is a dissipator describing the quality factor 
of the cavity phenomenologically as a photon amplitude 
damping channel [55]: 



-n)V[a] +nV[A] 



(2.9) 



where X'[^](p) = 2ApA^ - {A^A,p}, n = tr[Aapeq] 
characterizes the temperature of the bath, and k is the 
cavity dissipation rate (oc l/Q). The expectation value 
of the number operator at equilibrium is related to the 
temperature, T, of the bath by 



-1 

1) 



UJc 
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(2.10) 

where fcg is the Boltzmann constant. 

The reduced dynamics of the spin-ensemble in the in- 
teraction frame of the dissipator ( |2.9[ ) is given to 2nd 
order by the TCL master equation [26]: 



d , , 



dr tir 



Ci{t)e 



^lit-T)Ps{t)<g) P^g 

(2.11 



where Psit) = tvc[p{t)] is the reduced state of the spin- 
ensemble and Peg is the equilibrium state of the cavity. 
Under the condition that k 3> g^/Ws, the master equa- 
tion (2.11) reduces to 



d 
dt 



Psit) 



4 



t-to 



dre 



'^/^(^cos{AT)VsPs{t) 



-sm{AT)Csp{t) 



(2.12) 
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where 

Vs = {l + n)V[JJ\+nV[J+] (2.13) 

CsP = -i[Hs,p\ (2.14) 

H, ^ {l + n)J+J_-nJ_J+ (2.15) 

are the effective dissipator and Hamiltonian acting on the 
spin ensemble due to coupUng with the cavity. 

Under the assumption that k, ^ gV^s we may take 
the upper hmit of the integral in (2.11) to infinity to 



obtain the Markovian master equation for the driven spin 
ensemble: 



where 



dr 



Ps{t) ^ \ ilsCs + ^VA Ps{t) 



4A2' 



(2.16) 



(2.17) 



Here Jl^ is the frequency the effective Hamiltonian, and 
T s is the effective dissipation rate of the spin-system. 

If we consider the evolution of a spin state which is 
diagonal in the Dicke basis, p{t) — J2m ^rn{t)Pm, where 
Pm{t) — (^7 w|/9(t)| J, m) is the probability of finding the 
system in the Dicke state pm = |J, r7i)(J, m| at time i, 
the master equation (2.16) reduces to a rate equation for 



the populations of the Dicke states: 



d_ 
dt 



Pm{t) — T s(^Am+lPm+l{t) + BmPm{t) 

+ C„,-lP,n-l{t)) (2.18) 



where 



Am = (1 + n)[j(J + 1) - m(m- 1)] (2.19) 
Cm = n[j(J+l)-m(m + l)] (2.20) 

Bm = -{Am + Cm) (2.21) 



Defining P{t) = {P-j{t),. 
trix differential equation 



,Pj{t)), we obtain a ma- 



dt 



P{t) - TsMP{t), 



where M is the tridiagonal matrix 



M 



(B^j 

'C_.7 A_.7+2 

C_,7+i S-,7+2 



(2.22) 







V 



BjI 

(2.23) 

For a given initial state specified by Dicke-state popu- 



lations P(0), Eqn (2.22) has the solution 
P{t) = exp (ir, Af)P(O). 



(2.24) 



The equilibrium state of the driven spin-ensemble 
satisfies M ■ P{oo) = 0, and is given by p^ 



Y.'L=-j P'n{oo)pm, where 



-'eq 



P,„(oo) = 



=-./+m 



(1 + n) 



(1+fl) 

2,7+1 _ ^2.7+1 



(2.25) 



It corresponds to the non-thermal state 



exp 



T 



In 



1 



By comparing with Eqn (2.10) 



(2.26) 

we have that the final 
effective temperature reached by the spin ensemble in 
the interaction frame of the Rabi drive is T^q = ^T'c, 
where Tj, is the temperature of the cavity. The total spin 
expectation value for the equilibrium state of the spin- 
ensemble is 



-J 



(2 J -(- l)n 



2.J+1 



(1 + n)2.-'+i -n 



(2.27) 



In the limit of Ng ^ n, we have that the ground state 
population at equilibrium is given by P-j « 1/(1 -|- n) 
and the final expectation value is approximately {Jx)eq ~ 
— J + n. Thus, the final spin polarization will be roughly 
equivalent to the thermal cavity polarization. 

We note that if the detuning Slo were negative, match- 
ing flfi = 2Sui would result in the H^q^ term being dom- 



inant, leading to a master equation (2.16) with the oper- 
ators J_ and J+ interchanged, the dynamics of which 
would drive the spin ensemble towards the (Jx) = J 
state. Thus, the detuning must be larger than the cav- 
ity linewidth to prevent competition between the H^cir 
and ff+iifi terms, which would drive the spin system to 
a high entropy thermally mixed state. 



III. SIMULATIONS 



The tridiagonal nature of the rate matrix ( 2.23 ) allows 



Eqn (2.24) to be efficiently simulated for large numbers 
of spins. For simplicity we will consider the ideal case 
where the cavity is cooled to its ground state {n = 0), 
and the spin-ensemble is taken to be maximally mixed 
(PmiO) = 1/(2J + 1) for m = - J, . . . , J). 

The simulated expectation value of {Jx{t)} for ensem- 
ble sizes ranging from Ng = 10'^ to Ng = 10^ is shown 
in Fig. [T| normalized by — J to obtain a maximum value 
of 1. At a value of —{Jx{t))/J — 1 the total angular 
momentum subspace of the spin ensemble is completely 
polarized to the Jx ground eigenstate |J, — J). 

The expectation value {Jx{t)) may be fitted to an expo- 
nential to derive an effective cooling time-constant, Ti^off, 
analogous to the thermal spin-lattice relaxation time, Ti. 
A fit to a model given by 



J 



1 — exp 



t 



cff 



(3.1) 



4 




FIG. 1: Simulated evolution of the normalized expectation 
value of —{Jx{t))/J for a cavity-cooled spin ensemble con- 
taining A'^s spins. Time is in units of the effective dissipation 
rate Fs for the spin-ensemble given in Eqn ( pTf] ). 



FIG. 2: Effective cooling time-constant, Ti.off, of the spin 
ensemble as a function of the equilibrium excitation number of 
the cavity, n, for Ng = 10, 10^, 10"^, 10'* spins in the ensemble. 



yields the parameters Ti^off = XN'^/Ts with A = 2.0406 
and 7 = —0.9981. An approximate expression for the 
cooling time-constant as a function of the number of spins 
iV, is then 



Tl.eff(iV,) 



2{k^ + 4A2) 



(3.2) 



Eqn |3.2| shows that the cooling efhciency is maximized 
when the Rabi drive strength is matched to the cavity 
detuning (A = 0). In this case the cooling rate and time- 



g /k and Ti^, 



2K/g^N,, 



constant simplify to Tg 
respectively. 

If we assume a cavity cooled to its ground state with 
Q = 10^ (K/27r = 1 MHz) and a spin-cavity coupling of 
3/27r = 1 Hz j27j, the range of validity of the Markovian 
master equation is Ng ^ / 9^ — 10^^ and an ensemble 
containing roughly 10^^ electron spins may be completely 
polarized with an effective Ti of 3.18 /iS. This polariza- 
tion time is significantly shorter than the thermal Ti for 
low-temperature spin ensembles, which normally range 
from seconds to days [5]. 

In the case where the cavity is thermally occupied, the 
final spin polarization is roughly equal to the thermal 
cavity polarization, Eqn (2.261, and the fitted effective 



cooling constant of the spin ensemble as a function of 
the thermal cavity occupation number is shown in Fig [2] 
We find that for cavity temperatures corresponding to 
n < \/Ns the effective cooling constant Ti^^s is approx- 
imately equal to the zero temperature value. However 
when n > y/N^ the effective cooling constant is reduced, 
and approaches a value of Ti^off ~ (2fi)~^. 



IV. CONCLUSION 

The ability to reduce the effective Ti time of a spin 
ensemble by simply applying a detuned microwave drive 



provides an important tool for spin-based quantum in- 
formation processing (for example 28 - 30 and references 
therein), where fast gated reset operations provide con- 
tinual access to pure ancilla states for use in quantum 
error correction protocols. The ability to efficiently po- 
larize spin samples should also provide more reliable mea- 
surements in spectroscopy by permitting fast signal av- 
eraging. 

Several assumptions were made in the presented the- 
oretical model for cavity cooling of a spin ensemble. 
Firstly, we have assumed that the spin ensemble has no 
dipolar coupling between spins. This assumption is valid 
for magnetically dilute samples. The case of weak dipo- 
lar coupling between spins needs to be considered in more 
detail, but is not expected to fundamentally change the 
results. Secondly, we have neglected the effects of ther- 
mal relaxation of the spin system. As the cooling effect 
of the cavity on the spin system relies on a coherent spin- 
cavity information exchange, the relaxation time of the 
spin system in the frame of the Rabi drive - commonly 
referred to as Ti^p - must be significantly longer than the 
inverse cavity dissipation rate 1/k. Thirdly, we have as- 
sumed that the spin-cavity coupling, g, is spatially invari- 
ant and that the spin-cavity coupling and Rabi drive are 
spatially homogeneous across the spin-ensemble. Spatial 
inhomogeneities would lead to an incoherent distribution 
of cooling rates over the sample that obscures, but does 
not prevent, the cooling effect. 

Finally, the derivation of the Markovian master equa- 



tion ( 2.16 1 assumes that no correlations between the cav- 



ity and spin system accrue during the cooling process, 
such that there is no back action of the cavity dynam- 
ics on the spin system. This condition is enforced when 
the cavity dissipation rate, k, exceeds the rate of coherent 
spin-cavity exchange in the lowest excitation manifold by 
at least an order of magnitude - i.e. k > lOg^/N^ (see 
supplementary material). In this Markovian limit, the 
rate at which spin photons are added to the cavity is sig- 
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nificantly less than the rate at which thermal photons are 
added, meaning the cooling power of the fridge necessary 
to maintain the thermal cavity termperature is sufficient 
to dissipate the spin photons without raising the aver- 
age occupation number of the cavity. From eqn. (3.2 1 we 



see that, in principle, the cooling efficiency could be im- 
proved by adding more spins to make k closer to g\/Ns, 
but in this regime the cooling power of the fridge is no 
longer sufficient to prevent back action from the cavity 
and non-Markovian effects significantly lower the cooling 
rate (see supplementary material). 

As a final comment, the connection between standard 
techniques in quantum optics and magnetic resonance 



used in this work has been noted previously (for example, 
|31H36j ). but remains relatively unexplored. The many 
recent examples of coupling spin ensembles to high-Q 
resonators, microwave oscillators, and superconducting 
qubits (for example [571 and references therein) 

motivates more investigation into the relationship be- 
tween the two fields. 
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V. SUPPLEMENTARY MATERIAL 



In this supplementary material we provide additional details on the effective cooling rate for a thermally excited 
cavity and the validity of the Markovian assumption used to derive the cooling rates derived in the main text. 



A. Thermally excited cavity 



(3.2) 



In the Markov master equation simulations used to derive the effective cooling constant, Ti^off eqn. 
function of ensemble size, Ng, we assumed that the cavity was initially cooled to its ground state. If we fit our 
simulations for a non-zero average cavity occupation number n — (^q^i/^bTc _ ^ where oJc and Tc are the resonant 
frequency and equilibrium temperature of the cavity, the effective cooling time constant of the spin ensemble was 
found to obey 



l.off ' 



2 

J_ 

2n 



for n < y/Ng 
for n > iV, 



(5.1) 



This effect appears to originate from the fact that the final spin system polarization will be equal to the cavity 
polarization. Thus, cooling to a spin temperature that is not fully polarized requires removing fewer photons from 
the spin system. Given that the spin dissipation rate is independent of the cavity temperature, it takes less time to 
drive the spins to a state that is not fully polarized. An example of the simulated normalized spin expectation value 
~{Jx)/ J as a function of temperature is shown in Fig|3] Here we are considering a cavity with resonant frequency of 
Wc/Stt = 10 GHz for the Tavis-Cummings Hamiltonian in Eqn (2.3 1. 
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FIG. 3: Normalized spin expectation value —{Jx)/J of the spin ensemble as a function of time for various equilibrium temper- 
atures of the cavity. We consider the case of A'^^ =10, 100, 1000 and 10,000 spins in the ensemble, and a cavity with resonant 
frequency ujc ~ 10 GHz 



B. Validity of Markov Approximation 



The validity of the cooling rates derived by the Markovian master equation, eqn. (2.16), depends on the validity of 
the Markov approximation used to derive the master equation. As stated in the main text, enforcing the condition 
that the cavity dissipation rate, k, exceeds the rate of coherent exchange between the spins and cavity implies that 
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there will be no back action of the cavity dynamics on the spin system. We stated this condition as k > 10g-\/]V, such 
that any photons transferred to the cavity are immediately dissipated in the cavity before they return to the spin 
system. 
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FIG. 4: Energy level diagram of the joint spin-cavity system with coherent transitions denoted by a solid line and cavity 
dissipation rates denoted by a curved line. States are labelled as \ — + m,n) , where m is the number of spin excitations and 
n is the number of cavity excitations. For the cooling dynamics to appear Markovian states of high cavity excitation number 
should not be significantly populated on a coarse-grained time scale. 



More concretely, from the spin-cavity energy level diagram shown in fig. [4j the rate of transfer between states 
\ — + m,n) and \ — + ni — l,n + 1) is given by gy^(77ir+T)(iV^^^m), where N is the number of spins in the 
ensemble. At the same time, the cavity dissipator of strength ny'n is acting to drive the spin-cavity system to the 
state \ — Jx + m — l,n). To satisfy the Markov condition, we require the cooling dynamics to always drive the spin- 
cavity system toward states of low excitation number (bottom left of diagram) , without significantly populating states 
of high excitation number (top right of diagram). This will occur if the maximum rates for coherent transfer and 
cavity dissipation obey the following relationship: 

ffV(y + l)(A^-y)««Vy (5-2) 
Assuming that TV ^ 1, such that ^ + 1 « ^, this condition simplifies to 

K > O.mg^/N. (5.3) 

As shown in fig. [s] when n = 0.5gy/N a full non-Markovian simulation of the cooling procedure yields dynamics 
that are much richer than predicted by the Markovian model used in the main text. In particular, coherent transfer 
of spin photons deposited in the cavity back to the spin system are seen as oscillations in the expectation value of 
Jx- These memory effects reduce the cooling efficiency such that the cooling rate is initially very fast when the cavity 
occupation is low, then slows down significantly as higher excitations of the cavity are transferred back to the spin 
system. As k becomes larger than \/2, the oscillations are damped out, but the Markovian master equation still does 
not fully agree with the full non-Markovian simulation. When k = 10g\/N, the oscillations are critically damped and 
the Markovian master equation captures the full cooling dynamics. Thus, if the cavity dissipation rate, k, exceeds 
the rate of coherent spin-cavity exchange in the single excitation manifold by at least an order of magnitude - i.e. 
1^ ^ Wgy/N - then the Markovian master equation is valid. 
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FIG. 5: Comparison of the cooling dynamics by the Markovian master equation (pink dotted curve) and a full spin-ensemble 
cavity simulation (blue dashed curve). The normalized expectation value of —{Jx{t))/J is plotted for N — 5 spins with 
K — 0.5, 1,5, 10(?%/iVs and cavity temperature T = OK. When k > IQg^fN the Markovian master equation calculation agrees 
very well with the full simulation. Also, as predicted by eqn. (3.2 1, the cooling rate increases for larger until the point where 
non-Markovian effects take over. 



